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ABSTRACT: 

In this paper, the first microscopic approach to the Brownian motion is developed in 
the case where the mass density of the suspending bath is of the same order of magnitude as 
that of the Brownian (B) particle. Starting from an extended Boltzmann equation, which 
describes correctly the interaction with the fluid, we derive systematicaly via the multiple 
time-scale analysis a reduced equation controlling the thermalization of the B particle, i.e., 
the relaxation towards the Maxwell distribution in velocity space. In contradistinction to 
the Fokker-Planck equation, the derived new evolution equation is non-local both in time 
and in velocity space, owing to correlated recollision events between the fluid and particle 
B. In the long-time limit, it describes a non-markovian generalized Ornstein-Uhlenbeck 
process. However, in spite of this complex dynamical behaviour, the Stokes-Einstein law 
relating the friction and diffusion coefficients is shown to remain valid. A microscopic 
expression for the friction coefficient is derived, which acquires the form of the Stokes law 
in the limit where the mean-free path in the gas is small compared to the radius of particle 
B. 
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1. Introduction 

This paper is concerned with a microscopic theory of the Brownian motion performed 
by a massive particle suspended in a gas of much lighter particles. Because of the mass 
and length-scale difference between the two components, one expects to eliminate the fluid 
variables from the description of the system and obtain a closed equation for the B particle 
[1,2,3]. Traditionally, this task is done by assuming that the dynamical properties of the 
Brownian particle (B) evolve on a time-scale much longer than the characteristic time-scale 
of the fluid. This approach leads to the well-known Fokker-Planck equation, which governs 
the time evolution of the distribution function /(R, V;£) of the position R and velocity 
V of particle B : 

M denotes the mass. In this equation, the fluid enters only through the friction coefficient 
C and the temperature T. 

An equivalent description can be obtained starting from the stochastic Langevin equa- 
tion 

dV 

M— = -MCV(t) + F(t) (2) 

where V(£) is the velocity of B and F(t) is the fluctuating part of the force exerted by the 
fluid, assumed to have a white spectrum. This equation leads to an exponential decay of 
the velocity autocorrelation function with a relaxation time given by C 1 [1,2] . 

On the other hand, some attempts at a fully microscopic description of the dynam- 
ics have been made [4,5,6,7,8,9]. It has indeed been possible to get rid of the stochastic 
assumption and obtain the Fokker-Planck equation by a systematic expansion of the dy- 
namics of the complete system, bath + B particle, in powers of the square root of the mass 
ratio m/M, where m is the fluid particle mass. Such derivations lead in particular to a 
microscopic formula for the friction coefficient in terms of the autocorrelation function of 
the force acting on the B particle. 

However, as pointed out by many authors [10,11,12,13], but already by Lorentz in 1911 
[14], the description based on eqs. (1) and (2) can only be valid if the ratio p/ps of the 
mass density of the fluid, p, and that of the B particle, p^, is very small. Indeed, as 
noticed above, the velocity of particle B relaxes on a time-scale Ty ~ C 1 . If Stokes' law 
is assumed, this leads to ry ~ M/rjH, where r\ is the viscosity of the fluid, and £ the B 
particle diameter. On the other hand, a typical hydrodynamic time of the fluid is of the 
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order Tf ~ T, 2 /(n/p), which is the time for a shear flow to propagate over the distance E. 
These rough estimates give Tf/ry ~ p/pb, so that the assumption of a wide time-scale 
separation is only justified if this ratio is small. This condition is unfortunately far from 
the experimental situation, where the mass density ratio p/pb is taken rather close to unity 
to avoid sedimentation of particle B. In this case, the fluid dynamics must contain slowly 
decaying modes, which relax on the same time scale as the velocity of the B particle. In 
other words, non-Markovian effects are expected and the validity of the Fokker-Planck or 
Langevin equations becomes doubtful 1 [12]. 

Several attempts have been made to overcome these difficulties and determine the dy- 
namical evolution in the general case, that is for any (finite) p/ps [10,11,15]. The main 
idea underlying all these works is that the slowly decaying fluid modes result from the 
momentum conservation law for the fluid particles. A correct description should therefore 
treat both fluid and B particle variables on the same level. This can be done for exam- 
ple by using fluctuating hydrodynamics for the fluid motion, with appropriate boundary 
conditions on the surface of the suspended B particle [10]. These approaches lead to a 
non-Markovian Langevin equation, involving the time-dependent friction coefficient : 

r/V r l 

M— = -M J dr C(t - r) V(t) + F(t) (3) 

We present here the first systematic, microscopic approach which treats the B particle 
and the fluid at the same level, clarifying and providing the fundamental basis for such 
phenomenological descriptions. The problem is reconsidered within the kinetic theory of 
gases. The system consists of a single, large, heavy hard sphere (mass M, diameter E) 
suspended in a fluid of small, light hard spheres (mass m, diameter a). The dynamics of 
the system will be assumed to be governed by an extended Boltzmann equation, which 
correctly takes the gas-B particle collisions into account. The same type of equation was 
used by van Beijeren and Dorfman to develop the kinetic theory of hydrodynamic flows 
and in particular to clarify the dynamical origin of the hydrodynamic Stokes law for the 
friction coefficient [16]. The extended Boltzmann equation is expected to be correct in the 
Grad limit, defined for the gas with number density nbyn-^oo,cr^0, the mean free 
path I = (na 2 )^ 1 = constant [17]. We choose the diameter E of the B particle to fix the 

Let us note however that the Smoluchowski equation, describing the spatial evolution of 
the B particle is not affected by these arguments, since the position of the B particle relaxes 
on a time-scale much longer than the velocity (or fluid) relaxation times. 
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length scale. The condition of a constant mean free path thus imposes 

E o 

— = no E = const (4) 

Moreover, and unlike all previous work, we require from the begining the condition of 
essentially equal mass density of the gas and the B particle. In other words, the ratio 
p/PB wu l be kept finite even in the small mass ratio limit, which amounts to assume 

M 

^2 = const x mn (5) 
Finally, we define the small parameter corresponding to the Brownian limit 

e =(m) <<1 (6 » 

Introducing the diameter ratio k = cr/E, we characterize the asymptotic regime (4)-(6) 
by the limit : 

7iE 3 ~e~ 2 ( 7 ) 
To simplify the calculations and without loss of generality, we will assume in the following 
that e = k and nE 3 = e~ 2 . This can be looked upon as redefining the parameters of the 
system (E, M, n, ...), to set the two constants introduced in eqs. (4) and (5) equal to one. 
The correct scaling will however be recovered in the final results. 

Starting from the extended Boltzmann equation, we will perform a systematic expan- 
sion of the dynamical properties of the system in powers of e, in order to obtain a closed 
equation for the distribution function of particle B. 

The following results will be obtained : 

(i) On the shortest time-scale, t ~ ry ~ Tf, the B particle "does not move" (in position 
space), whereas its velocity distribution relaxes in a thermalization process. 

(ii) The equation (65) characterizing the relaxation of the velocity distribution is not 
a Fokker-Planck equation. It exhibits a complex non-markovian character, which stems 
from building up the dynamical friction force by recollision events between the gas and the 
B particle. 

(iii) We will be able to express the autocorrelation function (ACF) of the velocity of 
the B particle, in terms of the time-dependent friction coefficient. A microscopic formula 
for this coefficient follows from our analysis. The final result justifies the phenomenological 
expression for the velocity ACF, based on fluctuating hydrodynamics. 
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(iv) In spite of non-markovian effects, the Stokes-Einstein relation, expressing the 
diffusion constant of the B particle in terms of the friction coefficient, is explicitly derived. 



2. Kinetic Equations 

In the study of the dynamics of the B particle immersed in a bath of N small spheres, 
the short-hand notation will be used : 

S=(R,V); i = (ri,Vi), z = l,2,...,iV (8) 

for the positions and velocities of the B particle and the N gas particles. 

The temperature T of the gas introduces characteristic measures, ^/k B T/M and 
^k B T/m, of the velocities of the B and fluid particles, respectively. On the other hand, 
the characteristic length scale in the system is given by the diameter E of the B particle. 
In order to study the asymptotic regime (7), it is convenient to use dimensionless variables 
defined by 

R = EX, rj = Exj 

(9) 

The shortest time-scale in the system is fixed by the gas-B particle collision frequency, and 
we introduce accordingly a dimensionless time variable r by 



f ^V5 r (10) 

Our aim is to determine the dynamical evolution of the B particle distribution function 
f B (B; t). To this end, we need to introduce the joint distribution, f B i(B, 1; t), representing 
at time t the number density of fluid particles in state 1, when at the same time particle B 
is in state B. Finally, we define the conditional distribution function of the gas, f\\ B , by : 

f Bl {B,l-t) = f B {B-t) h\ B (l\B-t) (11) 

Under the dimensional scaling (9), the dimensionless distribution functions, F B and 
i*l|_g, are defined by the relations 

F B (B,r) = E 3 (k B T/Mf/ 2 f B (B;t) (12) 
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and 

F MB (1\B; r) = Y?{k B T/mfl 2 f l]B (l\B-t) (13) 

However, whereas F B (B;r) is essentially of order unity, F^ B does diverge as e — > 0. 
Indeed, far from the suspended sphere B, the conditional distribution function f^ B is of 
the order of the density n of the gas, which implies (according to condition (7)) : 

F^ B ~ n£ 3 ~ e~ 2 (14) 

Having in view a perturbation expansion, we thus introduce rescaled distribution functions 
F^ b (1\B;t) and F B i(B,1;t) defined as 

F MB (l\B;r) = e 2 F llB (l\B;r) 
F b1 (B,1;t) = 6 2 F b (B;t)F mb (1\B;t) 
which are of order unity when the small parameter e goes to zero. 

The scaled distribution functions evolve in time according to the coupled equations : 

^ + e 3 U- A]f b ( j B;t)= I dlT e _{B,l)F B1 {B,l-T\ (16) 



i h e J U • — + e z ui ■ 



T e _(B,l)^}F B1 (B,l;r) 



(17) 



Xdr ax V dxi 

rf2 {r (S, 2) + e 2 T_(l, 2)} F B12 (B, 1, 2; r) 

where the factorization of the three-particle distribution 

F B12 (B, \,2;t)~F B {B; t)F l]B (l\B- t)F^ B {2\B- 1) (18) 

is expected to hold in the Grad limit (see the discussion leading to (4)). 

Equation (16) expresses the fact that the B particle state changes in time through 
free-motion (left-hand side (l.h.s) of the equation) and through collisions with gas particles 
(right-hand side (r.h.s) of the equation). The second equation (17) states that the two- 
particle distribution, F B \, evolves in time, owing to free motion and collisions between 
the pair (B, 1) (l.h.s of the equation), and also owing to collisions of this pair with the 
surrounding gas particles (r.h.s of the equation). 

These evolution equations involve two hard-sphere collision operators, T < L(S, 1) and 
T_(l,2). The first of them characterises the effect of binary collisions between the B 
particle and a gas particle. It reads : 

T e _(B, 1) J da [(cU - ui) • a] 6 [(eU - u x ) ■ a] 

x {5(X-^-xi)6|(S,l)-5(X + ^-xi)} 
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where £ = 1/2 is the dimensionless radius of the B particle; the operator 6| transforms the 
pre-collisional velocities into post-collisional velocities : 

[6|(s, i) x ] (u, Ul ) = x (y- ^2 t( eU - u i) • *] ^ u i + rh 1 [(eU " Ul) ' a] a ) 

(20) 

where a is the unit vector along the axis joining the centres of the two spheres at contact 
and x is an Y function of velocities U and ui. 

On the other hand, T_(l,2) characterizes the effect of binary collisions between the 
fluid particles. It is given by : 

T_(l,2) = J da(u 12 -a)0(u 12 -a)5(x 12 ){^(l,2)-l} (21) 

where ui 2 = ui — u 2 , xi 2 = xi — x 2 , and acts as &| in eq. (20), with e = 1 
(collision between equal mass particles) . Note that the dimensionless cross section of the 
gas, (cr/E) 2 = e 2 (see eq. (7)), has been extracted from the expression of T_ and displayed 
as a prefactor in eq. (17). 

In eqs. (19) and (21), the fluid part of the collisional transfer is not present since we 
consider the Grad limit : at encounters, the positions of the particles coincide [18]. 

To solve the system of integro-differential equations (16)-(17), we have to define the 
initial state of the system. As in our previous work [8,9], we shall assume that the fluid is 
initially in conditional equilibrium in the presence of particle B : 

F e m (B, 1; r = 0) = F%{B- r = 0)F e ?(l|X) (22) 

with 

F e «(l|X) = e(|X-x 1 |-O0(ui) (23) 

where 

<P(u) = (27r)- 3 / 2 exp(-w 2 /2) (24) 

is the Maxwell distribution and © is the Heavyside step function. 

Before closing this section, let us recall some results concerning the e expansion of 
the collision operator T e _(B,l), which will be used in the course of the analysis. The 
calculations leading to them , though technical, are quite straightforward and have been 
presented in our previous paper (ref. [8] quoted later as I). 

The B-gas collision operator can be formally expanded in powers of e = (m/M) 1 / 2 as 

T e _(B,l) = T {0) (B,1) + eT^Bd) + e 2 T {2) (B,l) + ... (25) 
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The zeroth order term T^(B,1) characterizes collisions of the gas particles with the 
immobile B particle, acting as an external field. One finds 

T^(B, 1) =i 2 [ da (-ui • *)0(-ui • a) 

r >, (26) 

x U(X-t&-x 1 )bf\B,l)-6(X + &-xi)j 

where represents the change of velocity of a gas particle, undergoing a specular collision 
with the B particle fixed at point X (see eq. (20), with e = 0). 

The first order correction has a much more complicated form (see Appendix of paper I 
for the full expression), but we will only need the action of T^(B, 1) on a state in which 
the gas is in conditional equilibium (23). In this case, one finds the following formula : 

T { }\b, l)F(B)F e ?(l|X) = F(B)V • ^-F e <?(l|X) - .F_(l) • + Ij) F(B)F e ?(l|X) 

(27) 

where the notation T^{V) has been introduced to denote a dimensionless microscopic 
"force" 

We shall also need some formulas involving integration over the fluid degrees of freedom. 
First, one can verify that the integrated zeroth-order term identically vanishes : 

J dlT { °\B,l)F B1 (B,l) = 0. (29) 

The first-order term is given by 

f dlT W (B,l)F B1 (B,l) = 

J , d (30) 

= -J dlF+(l)-— F B1 (B,1) 

with ^+(1) defined in eq. (28). 

The general expression for the integrated second order term is more involved (see eq. 
(31) of paper I for the complete formula), but when the fluid is in conditional equilibrium, 
it reduces to 

J dl T {2) {B,l)F B {B)Fl q {l\X) 
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3. The Multiple Time Scale Analysis 

The system of coupled equations (16), (17) will be studied in the Brownian limit 
where e C 1. As in our previous work, this will be achieved systematically by using a 
multiple time-scale analysis, which leads to a uniformly convergent e-expansion, avoiding 
secular divergences as time goes to infinity. We just recall here the spirit of the method. 
In the e — > limit, different time-scales separate out in the system. Accordingly, we 
replace the distribution functions Fg, Fq\ by auxiliary functions, Fg(B; tq, ti, T2, . . .), 
Fg 1 (S, 1; to, 7~i, T2, . . .) which now depend on many time arguments. The time derivative 
is replaced accordingly by the operator 

+ e — + e 2 — + . . . (32) 



8tq dr\ 8t2 
The auxiliary functions are also expanded in powers of e : 

The expansions (32), (33) are substituted into the evolution equations (16) and (17), 
supplemented by (18), and terms of the same order in e are identified. The determination 
of successive corrections Fg\ F^ (k — 0,1,2,...) is then achieved by combining the 
chosen initial condition with the requirement that the expansion in e be uniform with 
respect to time, which amounts to eliminating secular divergences. The physically relevant 
solution of eqs. (16) and (17) is then obtained by restricting the originally independent 
multiple time variables tq, t\, T2 to the "physical line" : 

tq = t, , ri = er, ,T 2 = e t, ... (34) 

on which the operator (32) reduces back to -g^. The dependence of the distribution 
functions on the variable tj essentially defines the dynamical evolution on the time scale 
r — e ~K j = 0, 1, 2, ... 

Finally, we shall assume that the initial condition is entirely contained in the zero order 
terms : 

4 0) (£; to = o, n = o, r 2 = o, . . .) = f%(b- t = o) 

F%{ (B, 1; to = 0, n = 0, t 2 = 0, . . .) = F% X {B, 1; r = 0) 
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4. e-Expansion and Derivation of the Reduced Equation for the B Particle 

We first consider the zeroth-order terms in the coupled kinetic equations (16)-(17). 
Since at this order the integrated gas-B collision operator vanishes (see eq. (29)), we are 
left with 



—F^\B;t ,t 1 ,t 2 ,...) = 



dr 



—Fg}(B,l;r ,T 1 ,T 2 ,...) = 



(36) 



So both distributions, and F^, do not depend on the time variable tq. 
To first order in e, the kinetic equations reduce to 

The action of the collision operator T^) is given in eq. (30). Since F^ and F^ and the 
r.h.s. of eqs. (37) are independent of tq, we must impose 

' -F%\B;t ,t 1 ,t 2 ,...) = 



%° (38) 
— F^ L 1 ) (S,l;To,ri,T 2 ,...) = 

to eliminate secular divergences. On the other hand, one can verify that eq. (30) implies 

dl T^\b, l)F B (B)F e %l\X) = (39) 



The integrated operator applied to a state of the system in which the fluid is in 
conditional equilibrium, vanishes after integration over the fluid variables, whatever the 
state of particle B. Hence, the solution of eqs. (37) consistent with the initial condition 
(22) has the form : 

A F (0) (jB;Tl , T2 ,...) = 



(40) 

(1\B) = F e «(l]X) 

The evolution of the system takes thus place on a longer time-scale, corresponding to 
variable r 2 . The second order terms in the kinetic equations (16)-(17) yield the relations 

irf), iji), l F w_ 

8 ^ B , " * " , _ (41.a) 

= dl T ( l\b,1)F { ° ) F e %l\X) + dl T ( L ) (B,1)F { 2 ) 
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= ^2 T i2) (B,2)F { ® ) F eq (l\X) F eq (2\X) + J d2 T ( l\b,2)F { ^ ) F eq (l\X) F, 
+ Jd2 T_(1,2)F^ 0) F e9 (l|X) F e,? (2|X) 



(41.b) 

In obtaining these equations, the equality = F eq {l\X) and the relation (39) have 

been taken into account. 

Since F^ and F^ are independent of tq and t\, and and Fgj of to, we conclude 
from (41) that in order to eliminate secular divergences one has to impose 

(42) 



— Fg\B-T ,T 1 ,T 2 ,...) = —Fg}(B,l-T ,T 1 ,T 2 ,...)=0 



A F a) (jB;T1 , T25 ...) = JL F g}(B,i-,T h T 2 ,...) = A F w ( i| S;T1)T2 ,...) = 

Now, since -F e<? (l|X) is the conditional equilibrium of the gas in the presence of the 
fixed B particle, the following relations are satisfied : 

(ui • ^ - T^(B, 1)) F eq (l\X) = (43.a) 

J dl T_(l, 2)F e < ? (l|X)F e < ? (2|X) = (43.b) 

In view of these equalities, both evolution equations (41. a ) and (41. b ) are found to reduce 
to 



dr 2 
da ( 

coefficient by 



a - F p = c B ^ ■ (A + u) + / a !)<» ,44) 



where formula (31) for T_ (5,1) has been used and we defined the Boltzmann friction 

Cfl = e 2 |v^F (45) 



The first term on the r.h.s of (44) characterizes the instantaneous "static" friction force 
induced by collisions of the gas particles with particle B. By analogy with the hydrodynamic 
expression for the friction force (£x velocity), the role of the velocity is played here by 
(mJ ^) Fb\b). On the other hand, the physical meaning of the second term on the 
r.h.s. of (44) can be clarified by introducing a mean dynamic friction force F+ defined by 

T+{B-r 2 ) = Jdl F + (1)F^(1\B;t 2 ) (46) 
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This time-dependent force characterizes the "dynamical" part of the drag on the B particle, 
induced by the correlations building up inside the gas, in response to the motion of particle 
B. Using this definition, eq. (44) can be rewritten 

£4Vi • • ■) = £j ■ {<B (± + U) Ff(B; „, . . .) - 7 +{ B; r 2 )F^B; r 2 , . . .) 

(47) 

To close this equation, we have to evaluate the friction force (46). To this end we must 
consider the third order terms in the kinetic equations (16), (17). 

We first note that the already presented reasonning concerning the secular divergences 
(see eq. (42)), when applied to the third order equations, yields 

— 4 3) (S;t ,ti,t 2 ,...) = — Fg}(B,l;T ,T h T 2 ,...) = 



(48) 

The evolution equations to third order in e take thus the form 



S-F§\b; Ti ,t 2 ,...) = ^F§l(B,l-r h T 2 ,...) = ^-F 1 ( g(l|S;r 1 ,r 2 ,...) = 



= dl T a) (S, 1)^ + / dl T {2) (B, 1)F$ + dl T {3) (B, 1)F { ®I 
^{^4) + ^)^(l|X)} + ^)^(l|X) 

+ (ui • ^ - T?>(fl, 1)) F<°> F« + {U ■ A _ 5ll) (fl> 1} } F (0) F e, (1|x) 
= y d2 tL 1) (5,2)F^ 2 1 ) 2 + J d2 T ( L ) (B,2)F%} 2 + J d2 T < L ) (B,2)F^I 2 
+ Jd2 T_ (1,2)4^ 

(49.b) 

where the relation (43) has already been taken into account. The functions F^ = 
{FgF 1 | 5 }^) and F^ 2 = {FbF\\b^2\b)^ in the collision terms represent the fcth or- 
der contributions /c = 0, 1, 2. For instance, 

F Bi2( B i !' 2 ) = { F S F 2\b}^ ^ 

\X) 

(50) 



= F^ 2) F e9 (l|X)F e? (2|X) 



+ F e *(2|X) + F e *(l|X)F«} 

+ {f^ F e ?(2|X) + F e ^(l|X)F 2 (2 ] + F« F« } 
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Our concern is now to obtain a closed equation for the conditional distribution of the 
gas F^g(l\B). To this end, we have to eliminate the terms containing the T3- variable, 
which can be achieved by subtracting from (49. b ) the first equation (49. a ) multiplied by 
F e <?(l|X). One finds 

9 AV\ , 9 



dr 2 



«)^») + (u 1 -^-T?)( B ,l))^)^) = 

= T a) (S, 1)F^ F e ?(l|X) - i^ 0) U • ^F e %l\X) 

+ Jd2 T { 1\B, 2)4 0) F« (51) 

+ J d2 T ( l\b,2)F { ® ) Ffy F e ?(2|X) 

+ Jd2 T_(l,2)4 0) {f® F e ?(2|X) + F e *(l|X) F®} 

The last term of (51) introduces the (dimensionless) linearized Boltzmann operator A_g(l), 
defined as [18] : 

A B (1)»(1|B) = Jd2 F-(l|X)^(u 2 ) T_(l, 2) { (™) + (£g*>) } (52) 

where 0(u) is the Maxwell distribution. Note that A_g(l) acts only on the fluid variables. 
With the use of relations (27), (30) and (31), the equation (51) reduces to : 

= --Mi) ' (^j + u ) fiV)f e? (i|x) 



The notation 



£ / = u 1 .^--tL 0) (s,i)-a b (i) 



(53) 



(54) 



and 



F corr (B, l;r 2 ) = F^(B)F^ (1\B), (55) 



permits to cast (53) in a more transparent form 

{A + C f + £ B } F C0 -(S, 1; r 2 ) = -.F_(l) • (A + u) Ff(S; r 2 )F^(l|X) (56) 
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The force T ' j, appearing in the definition (54) of £5, can be interpreted as the total friction 
force acting on particle B. The system of equations (47), (56) is closed and characterizes 
the dynamical evolution of the distribution functions on the t 2 time-scale. 

The formal solution of eq. (56) can be written as 

F corr (B, 1; r 2 ) = - jH ds exp j- jH da' {C f + C B ) (*')} 



(Mj +U )^ 0)(jB;S)jPe,?(1|X) 



(57) 



However, as can be verified from the definitions (54) combined with relations (26), (52) 
and (46), Cf acts only on the fluid variables, while L B acts only on the velocity of the B 
particle. The two operators thus commute and the solution for F corr (B, 1) reads 



F C0 -(S, 1; 72) = - f 2 ds ^_(1; -(72 - a))F^{l\X) 
JO 



exp {- jH ds' C B {s>) } (^J + u) Fjp(B; s) 



(58) 



where J--{1; — r) denotes the force ^"-(1; 0) on particle B propagated by the intrinsic fluid 
dynamics in the presence of the B particle, fixed at point X (see the definition of 
(26)), backward in time to the instant —r. 

The dynamical part of the mean friction force, JF + (S), defined in (46), can now be 
obtained by multiplying (58) by ^+(1) and averaging over the fluid variables. This yields 

T + {B-t 2 )F b Q \b-t 2 ) = J dlF + (l)F corr (B,l;T 2 ) 

= J T J ds C d y n (r 2 - s) exp j - ^ ds' C B (s') } ( A + F b °\b- s) 

(59) 

where we have introduced C^ n (r), a force-force correlation function, defined by 

C dyn (r) = \ [ dl^-(l;-r).^ + (l;0)F^(l|X) 

I J (60) 

= -(^_(l;- r ).^ + (l;0)) MX) 

By using the notation (. . •) ( e g|X) ' we wou ld like to stress the fact that the average is taken 
over the equilibrium ensemble of the fluid, in the presence of particle B fixed at point X. 
The prefactor 1/3 stems from the isotropy of the system around the B particle. 

Moreover, eq. (59) can be seen as a self-consistency relation for the (time-dependent) 
mean force JF + , because the propagator C B , defined in (54), depends itself on the friction 
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force T+. Substituting this expression for T+ in eq. (47), one arrives at a reduced equation 
for the distribution of the B particle 



^ 0) (B;.)^.(A + U)^;.) 

jH ds C d y n (r 2 - S ) A • exp {- jH ds' £ B (s')} (A + u) ^(S; a) 



8t2 

+ 



where C B is defined in (54) and C dyn in (60). 

At this stage, we can return to the original variables R, V and t, using relations (9), 
(10), (12) and (13). For the sake of simplicity, we keep the same notations for ( B , C dyn , 
and for different forces and propagators involved in the evolution equations. Their full 
expressions are now given by 

CB = Tt{l) \< 2 ™ k BT) 1/2 (62) 

= SM^T^-fr"*) • f+(1; ° )> (* (63) 
where the formula for the microscopic forces T± now reads 

= (|) J da 2m(v 1 -a) 2 e(Tv 1 -a)a 6 (it- a - (64) 

In (63), the notation (eg|R) stands for an average over the equilibrium ensemble of the gas, 
in the presence of particle B, fixed at R. The dynamics is characterized by the propagator 
£s, given in the appendix together with C B (for the dimensionless form, see eq. (54)). 

Keeping in mind that we are only interested in the dynamics of the system on the time- 
scale characterized by the variable T2, i.e. r ~ e~ 2 (r is the dimensionless time variable 
defined in eq. (10)), we obtain the following equation for the distribution function f B of 
the B particle 

|/*(W B ^(v + ^A) /B(B;t) 

+ f ds C*»« - .) ± • exp {- f ds' C B (s>) } (V + tjf A) MB; s) 

_(65) 

This should be completed by the self-consistency equation for the dynamical part T+ of 
the friction force 

T+(B- t)f B {B- 1) = J* ds C d y n (t - s) exp |- jT* ds' C B (s') } (V + ^ ^ s) 

(66) 
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The system of equations (65)-(66) is closed and represents the dynamical evolution of 
the state of the B particle on the r ~ e~ 2 time-scale. This condition can be rewritten 
e 3 r <C 1, so that using eq. (10), we conclude that eq. (65) applies for times t verifying 



The spatial diffusion process has not yet started at this time scale, since time t is very short 
compared to the time needed to cover B particle radius with thermal velocity. However 
the fluid hydrodynamics is already at work. The time-scale (67) thus characterizes the 
relaxation of the velocity of the B particle, while its spatial state is not yet affected and 
will relax only on a longer time-scale. 

Clearly, the derived equation is not of a Fokker-Planck type. The main reason is that 
the friction force due to the bath builds up on the same time-scale as that characterizing 
the gas dynamics, which leads to memory effects in the relaxation of particle B. The friction 
force J- ' j has thus to be constructed in a self-consistent way and depends on the whole 
history of the Brownian motion. 

Equation (65) and its systematic derivation are the main results of the present work. 
5. Long-time Limit of the Reduced Equation 

In this section, we discuss briefly the long-time limit of the system (65)-(66). Note 
however that we stay in the time window r ~ e~ 2 , since it is only there that these equations 
do represent the evolution of the system. On longer time-scales (e.g. r ~ e~ 3 , r ~ e~ ), 
we expect the spatial relaxation to occur through the Smoluchowski equation (see the note 
placed in the introduction) [12]. 

In the long-time limit, the B particle velocity distribution will relax towards the sta- 
tionary solution of (65), that is towards the Maxwellian distribution. In the final stage, 
we can then write 



t < 




(67) 




(68) 



and the friction force J- j will accordingly decay to zero too. Then to first order in this 
quantity, we can put ~ in the expression (66) for JF + and obtain 
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The limiting form of the reduced equation for fs(B; t) then simplifies to : 




(70) 



where the time-dependent friction coefficient 



m = Cb s(t) + c d y 



(t) 



(71) 



has been introduced. S(t) is the Dirac distribution. 

This asymptotic form for the reduced equation of the B distribution calls for several 
comments. First, though much simpler than the complete equation (65) (together with 
(66)), it still exhibits a non-markovian nature. However, the non-locality of (65) in velocity 
space -hidden in the self-consistent definition of JF + in (66) - is now removed and the 
simplified form (70) only involves the Fokker-Planck operator ■ (v + ^ff-^y-)- In this 
sense, the velocity can be seen as following a "generalized Ornstein-Uhlenbeck process", 
still characterized by a transition probability sharply peaked around the mean value (i.e., 
only small jumps occur), but now keeping the memory of its whole previous history. 

In particular, whereas for the Ornstein-Uhlenbeck process the B velocity relaxes in an 
exponential way [1,3], the long-time behaviour of the friction coefficient ((t) modifies the 
nature of the thermalization process : we expect f^(B;t) to decay in an algebraic way 
towards the Maxwell distribution, since £(£) is known from hydrodynamic arguments to 
exhibit a t~^/2 long time tail [19]. This behaviour can be explicitly verified on the exact 
solution of eq. (70), which can be found in refs. [20,21]. 

6. Velocity Autocorrelation Function and the Stokes-Einstein Relation 

As we discussed above, the reduced equation (65) (together with (66)) characterizes 
the relaxation of the velocity of the B particle, while the spatial relaxation will occur on 
a longer time scale. As we shall see in this section, this time-scale separation allows one 
to compute explicitly the velocity ACF of the B particle. In spite of the non-local nature 
of the complete system (65)-(66), we will eventually recover the Stokes-Einstein relation 
between the diffusion and the friction coefficient [2]. 

Let us study the dynamical evolution of the mean velocity of the B particle, defined 

as : 




(72) 
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Since Vg(t) relaxes on the time scale r ~ e~ 2 , /g will be assumed to evolve according 
to equation (65). The evolution equation for Vg(f) can then be obtained by multiplying 
(65) by V and integrating over the velocity. This yields 



> = /^C 5 vA.(v + fA)/ B(B ;t) 

+ J* ds C d ^(t -s) J dVV A . exp j- jf* du £ B (u)} (v + f B (B; s) 

(73) 



The first term is calculated by an integration by parts, yielding 



dV vA.( V+ ^:^_ )/i}(B) = _v i} ( t ) (74) 



iV:(h; =,-,,,{ - H.£ B HUv + ^J-] /s ( ft „ ,75) 



The second term needs a more careful analysis. Let us introduce the notation 

M dV 
We thus have to compute 

/ vA. 7(V ;t| s ) = - J dV 7 (V;t|s) (76) 

where an integration by parts has been performed. According to (75), 7(V; t\s) solves the 
following initial value problem 

(! + £ B (t)) 7 (V;t| S ) = t>s 
l(V;t = s\s)=(v + ^^)f B (B;s) t= 77 



s 



d_ 
di 



Integrating eq. (77) over V, one obtains 

J dV 7 (V;£|s) + J dX C B (t) 1 (V;t\s) = (78) 

But according to eq. (A.l) in the appendix, the action of the operator C B involves multi- 
plication by the total friction force T j and then derivation with respect to V. The second 
term in (78) thus vanishes, since 

J dV £ fl (f) 7 (V;f| a ) = JdV A • {Tp(V;t\s)} = (79) 



and we are left with 



-| f dV 1 (V;t\s) = (80) 
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We thus eventually find 

J d\ 7 (V; t\s) = J dV 7 (V; t = s\s) 

= V B {s) 

Combining eqs. (73), (74), (76) and (81), we obtain a closed equation for of the form 

dt jo 

where ((t) = Cb^) + C dyn (t) is the time-dependent friction coefficient. 
Introducing the Laplace transform 

»+oo 



-V B (t) = -J*ds((t-s)V B (s) (82) 



r-too 

(z)= / dt exp(-z -t) V B (t), (83) 
JO 



one obtains the explicit solution for V B (z) as 

v bW = X« (84) 

where ((z) denotes the Laplace transform of the time-dependent friction coefficient. 

Equation (84) characterizes the relaxation of the mean velocity of the suspended sphere 
in a situation in which the B particle is initially out of equilibrium. But according to 
Onsager's principle of regression of fluctuations, this equation should as well describe the 
relaxation of a fluctuation of the velocity of the B particle at equilibrium [22] . This can be 
explicitly verified here. 

At equilibrium, the B particle velocity autocorrelation function is defined as 

(V(t) ■ V(0))eq = J dBdl . . . dN V(t) ■ V(0)p eq (B, 1 . . . N) (85) 

where p eq (B, 1 . . . N) is the canonical equilibrium probability density of the system, and 
V(t) denotes the velocity of the B particle propagated in time through the dynamics of 
the complete system. Formally, this relation can be rewritten 



(V(t) ■ V(0))eq = J dB V(0) • j dl . . . dN V(t)p eq (B, 1 . . . N) 
= JdB V(0) • (V(t)) n . e F eq (V(0)) 



(86) 
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where F eq (\) is the equilibrium distribution function of particle B alone, and (V(t)) n . e . 
denotes the time-dependent mean velocity of the B particle, averaged over the fluid vari- 
ables, for a given initial out-of-equilibrium state of velocity V(0). 

But this mean velocity will verify the evolution equation (82), derived previously by 
eliminating the fluid variables. This yields for the Laplace transform of the velocity ACF : 

ML 

(V(z) ■ V(0)) e , = -M- (87) 
z + ((z) 

The diffusion coefficient can be directly deduced from this relation, using : 

roc 

D S J o *(VW-V(0)U (8g) 
= (V(z = 0) ■ V(0)) e , 

This yields 

k B T 
D= M 



C(z = o) (89) 

M( f 

where, according to (71), the microscopic expression for the (integrated) friction coefficient 
Cf is given by 

/0 (90) 



r-\-oo 

Cf = ^ dr C(t) 

i r°° 



with given in (62), and the microscopic forces defined in (64). 

Equation (89) is precisely the Stokes-Einstein relation for the friction coefficient. 



7. Friction Coefficient and the Hydrodynamic Limit 

In this section, we would like to discuss the microscopic formula (90) obtained for the 
friction coefficient. Our aim is twofold. First, we would like to connect these results to 
the friction coefficient showing in the response of the fluid to an imposed motion of the B 
particle. And secondly, we would like to obtain, at least in some limits, explicit expressions 
for the friction coefficient as a function of the microscopic characteristics of the system (B 
particle diameter, transport coefficients of the fluid, . . .). This can be achieved by using the 
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results of van Beijeren and Dorfman (DVB), concerning the kinetic theory of hydro dynamic 
flows [16]. 

In accordance with (90) , the friction coefficient is the sum of two terms. The first term, 
£g, characterizes the static effect of instantaneous collisions between the B particle and 
the gas. The remaining part of the friction coefficient involves a time-integral of a force- 
force correlation function, reflecting the dynamical correlations induced by time-displaced 
collisions between the gas and the B particle. As we discussed in our previous work (ref. 
[23] quoted later as II), formula (90) is equivalent to the Kirkwood formula for smooth 
potentials, relating the friction coefficient to the time-integral of the force autocorrela- 
tion function. In the case of hard-sphere interaction however, the force is replaced by the 
momentum transferred to the B particle during instantaneous collisions and the separa- 
tion into the static term and the dynamical part then occurs (see II for further details). 
Moreover, it is interesting to note that the same microscopic expression for the friction 
coefficient was obtained in the limit considered in our previous work : e — > 0, all other 
parameters kept constant. As discussed in the introduction, this limit leads to the Fokker- 
Planck equation for the B particle, so that its evolution exhibits a Markovian nature. We 
thus see that the same properties of the fluid are involved in both limits, but they enter 
the dynamical laws at different levels. They are intrinsic hydrodynamic properties of the 
fluid. 

Let us consider the case in which the B particle has an imposed velocity U(£). This 
amounts to represent the B particle distribution by a Dirac 5-function in velocity space, 
centered on the mean velocity Vg(t) = U(t), with no thermal fluctuations around this 
value. The mean friction force T ' j acting on B can be then obtained from eq. (82), 
determining evolution of the mean velocity of the B particle 

7 fit) = - J^ds M((t - s) U(s), (91) 

which identifies £(£) as the time-dependent friction coefficient in the hydrodynamic sense. 
We would like to stress the fact that this result holds (as shown in section 6) although the 
reduced equation (65) is non-local in velocity space ! 

Let us now consider the dynamical part of the friction coefficient, which we defined in 
(71) and (90), as 

Cdyn(t) = ^^ f J dl ^rf S ^_(l;- S ).^ + (l;0)/ e ^(l|R) (92) 
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Our aim now is to connect this expression to the calculations of DVB, who have computed 
the drag on a macroscopic sphere, moving with constant velocity. The surrounding gas 
was assumed to obey an extended Boltzmann equation. 

To this end, we first introduce an auxiliary function / C (1|R), defined as 

rt 



/ c (l|R;t)= / dsF-(l;-s).V f%l\R) 

Jo t (93) 
= J ds exp(-s£ / )J !r _(f) • V f eq {l\R) 



where Cj is the fluid propagator (see eq. (A. 2) in the appendix). V is a unit vector intro- 
duced to establish the connections with the DVB analysis. According to (93), / c (l|R;t) 
verifies the following differential equation : 

(I + C f) /C(1|R; t] = ' ^ (94) 

Now if we define \P(l|R;t) through the relation 

/ c (l|R;t) = / e «(l|R) *(l|R;t) (95) 



the equation (94) can be rewritten 



(^ + £ / )*(l|R;t) = -T , (B,l) 



m 



vi • V (96) 



k B T 

where vi is the velocity of a gas particle, and the operator T'(B,1) is that defined in 
DVB's work : 

T'(B, 1) = (|) J d&i-V! ■ *)0(-vi • a)S(R - |a - n) [bf\l3, 1 ) - f] , (97) 

the operator changing the velocity vi of the gas particle to — 2(vi • a)a. 

Equation (96) is equivalent to the inhomogeneous Boltzmann equation (4.4) in the DVB 
paper. In their work, \I/(1|R;£) represents the dynamical correction to the distribution of 
the gas, induced by recollision events between the gas particles and the suspended sphere, 
which moves with velocity V. The authors have obtained approximate solutions of this 
equation in the limit where the mean- free path t is small compared to the radius S/2 of 
the suspended sphere. This solution has been obtained by decomposing the distribution 
functions into a hydrodynamic part and a boundary layer part. We refer to their paper for 
further details [16]. Their solution may be used to compute the dynamical part of Cdyn- 

22 



Indeed, with the aid of (92), (93) and (95), the dynamical contribution can be rewritten 



as 



Mk B T( dyn (t) = J dl F+(l;0)*(l\R;t) fl(l\ll) - V (98) 

where we used the isotropy of the system when introducing V. In DVB, the r.h.s. of (98) is 
interpreted as the dynamical part of the drag force exerted by the fluid on the suspended 
sphere, moving with velocity V. The established connection between both approaches 
reflects the linear-response theory, and can be physically interpreted through Onsager's 
principle of regression of fluctuations. 

In the limit of long times, we can then use the results (6.43) and (6.44) of DVB, to 
obtain to lowest order in i/T, 

Cdyn(t^00) = ^^-C B (99) 

where ( B is the Boltzmann friction coefficient, defined in (62), and r\ the shear viscosity of 
the gas, estimated within the Boltzmann approximation. Collecting both contributions C,q 
and Cdyn(t ~^ °°) > we then recover, as expected, the Stokes law for the friction coefficient 

poo 

c f = dt at) 

Jo (100) 
~ M 

In this limit, the diffusion coefficient acquires its Stokes-Einstein form 

2tt?7£ 



8. Conclusions 

In this paper, we considered the Brownian motion of a single heavy particle moving in 
a bath of light particles. Our study started with a microscopic description of the system, 
in terms of the coupled dynamical evolution equations for the distribution functions of 
the B particle and of the host fluid. The latter was assumed to evolve according to an 
extended Boltzmann equation, which correctly describes the effects of collisions between 
the gas particles and the suspended sphere. 
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Our aim was to derive a reduced equation for the B particle, by eliminating the gas 
degrees of freedom in the limit of small mass ratio e = (m/M) 1 / 2 . However, even in the 
e — > limit, we maintained the mass densities of both components to be of the same 
order of magnitude. The multiple time-scale analysis has been used to construct a uniform 
expansion in e. We derived in this way a new reduced equation (65) , governing the evolution 
of the velocity distribution of particle B. This equation turned out to be non-local in time 
and in velocity space. The corresponding memory terms result from building up of the 
friction force by the reaction of the suspending gas to the motion of B. However, in spite 
of its non-markovian character, this reduced equation allowed to compute the velocity 
autocorrelation function of the B particle (87). Moreover, we recovered explicitly the 
Stokes-Einstein law, relating the (spatial) diffusion coefficient to the friction coefficient 
(89). Finally, the derivation yielded a microscopic expression for the friction coefficient, 
reducing to the Stokes law in the hydrodynamic limit, where the mean-free path of the gas 
is small compared to the radius of the suspended sphere (100). 

Our work is not the first attempt to obtain the reduced equation for the distribution 
function of the B particle from a microscopic point of view [4,5,6,7,8,9]. However, all previ- 
ous approaches derive a (local) Fokker-Planck equation governing the dynamical evolution 
of the B particle by taking the m/M — > limit with all other parameters kept constant. 
Such a limit implies the asymptotically vanishing mass density ratio, while in the present 
work, we maintained this ratio at a constant value. This crucial difference is the source of 
the non-markovian character of the reduced equation for the velocity distribution of the 
B particle, since the velocity decays then on a hydrodynamic time-scale of the suspending 
bath. 

The results derived here suggest some directions for future work. First, the extended 
Boltzmann equation was taken as a starting point, which can be expected to be valid 
in the Grad limit. This conjecture could perhaps be verified explicitly by starting from 
the BBGKY hierarchy for hard-spheres [18], and applying directly to it the present limit 
defined in (7). In particular, we could in principle compute the e expansion of correlations 
induced in the bath by the presence of particle B. It would be then possible to check 
whether the results obtained in this work become indeed exact in the limit (7). 

On the other hand, it would be interesting to generalize the analysis to an arbitrary 
number of suspended particles, in order to obtain the generalized reduced equation for the 
N particle suspension. In this case, we expect the hydrodynamic interactions between the 
different suspended particles to be non-local in time too. 
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The reduced equation we obtained here characterizes the relaxation of the velocity dis- 
tribution of B, while the spatial distribution evolves on a much longer time-scale. Because 
of this time-scale separation, one then expects the spatial distribution to relax according 
to the local Smoluchowski equation. In principle, one could recover this equation from the 
same microscopic approach by pursuing the expansion up to the corresponding time-scale, 
t ~ E 2 /D, which corresponds to the T4 variable. Because of the complexity of the anal- 
ysis, we have performed up to now the corresponding calculations only for the case of an 
ideal (collisionless) suspending gas. Even in this case, in spite of the absence of hydrody- 
namic modes in the bath, memory effects are still present and they show once two or more 
Brownian particles are present. The work along these lines is being carried presently. 
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APPENDIX 



In this appendix, we give the full expressions for the propagators and Cf. 
With proper dimensions, can be written 

£B(0 = ^(c S (v + ^)-F + ( B ; t) ) (A,) 

where the mean force (.£?;£) is defined self-consistently by eqn. (66). On the other 
hand, the fluid propagator Cj reads : 

C f (l) = vi ~ - T^V, 1) - A B (1) (A.2) 

where A_g(l) is the linearized Boltzmann operator, and T^(B, 1) characterizes the effect 
of collisions between and the gas and the B particle, fixed at point R. The expression of 
A_g(l) is given by [18] : 

A B (1)¥(1|B) = f eq (l\n) Jd2 J da [(vx - v 2 ) • a] 6 [( Vl - v 2 ) • a] (5(n - r 2 )0(v 2 ) 

*«^-M(»ws?)} 

(A.3) 

where ^>(1\B) is a function of the fluid variables 1 = {i"i,vi}, and 0(v) is the Maxwell 
distribution 

The operator 6^.(1,2) acting on a function x(vi,v 2 ) replaces the velocities (vi,v 2 ) by 
their post-collisional values 

[hO-i 2 )x] (vi, v 2 ) = x (vi - ((vi - v 2 ) • a) <t, v 2 + ((vi - v 2 ) • a) a) (A.5) 
Finally, T^(B,1) can be written as 



^(S, !) = (§) / [-vi • a] [- Vl ■ a] 



x <; 5 ( r - |«r - bf\B, i) - s(k + y - n 



(A.6) 



where b^P\B, 1) changes the velocity vi of the gas particle into V]/ = vi — 2 [vi • a] a. 
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